Partial dynamical symmetry as a selection criterion for many-body interactions 
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We propose the use of partial dynamical symmetry (PDS) as a selection criterion for higher-order 
terms in situations when a prescribed symmetry is obeyed by some states and is strongly broken in 
others. The procedure is demonstrated in a first systematic classification of many-body interactions 
with SU(3) PDS that can improve the description of deformed nuclei. As an example, the triaxial 
features of the nucleus 156 Gd are analyzed. 
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Many-body forces play an important role in quantum 
many-body systems They appear either at a funda- 
mental level or as effective interactions which arise due to 
restriction of degrees of freedom and truncation of model 
spaces. A known example is the structure of light nuclei, 
where two-nucleon interactions are insufficient to achieve 
an accurate description and higher-order interactions be- 
tween the nucleons must be included Q. Given the dif- 
ficulty in constraining the nature of such higher-order 
terms from experiments, one is faced with the problem 
of their determination. One way, currently the subject of 
active research Q, is to determine them from chiral ef- 
fective field theory applied to quantum chromodynamics. 
This establishes a hierarchy of inter-nucleon interactions 
according to their order. In light-medium nuclei, these 
interactions serve as input for ab-initio methods {e.g., the 
no-core shell model (NCSM) 0|) to generate, by means 
of similarity transformations, A-body effective Hamilto- 
nians in computational tractable model spaces. 

The situation is more complex in heavy nuclei, where 
ab-initio methods are limited by the enormous increase 
in size of the model spaces required to accommodate cor- 
related collective motion of many nucleons. One possible 
approach to circumvent this problem, is to augment the 
NCSM method through a symplectic symmetry-adapted 
choice of basis [j| . A second approach is to employ energy 
density functionals and incorporate beyond mean- field ef- 
fects by mapping to collective Hamiltonians Q, e.g., the 
interacting boson model (IBM) [7]. In both approaches 
the Hilbert spaces are based on particular dynamical al- 
gebras which lead to a dramatic reduction of the basis 
dimension. Nevertheless, even with such simplification, 
the number of possible interactions in the effective Hamil- 
tonians grows rapidly with their order, and a selection 
criterion is called for. In this Rapid Communication, we 
suggest a method to select possible higher-order terms 
which is based on the idea of partial dynamical symme- 
try (PDS). 

The concept of PDS @ is a generalization of that of 
a dynamical symmetry (DS) H where the conditions of 
the latter (solvability of the complete spectrum, existence 
of exact quantum numbers for all eigenstates, and pre- 
determined structure of the eigenfunctions) are relaxed 



and apply to only part of the eigenstates and/or of the 
quantum numbers. PDSs have been identified in various 
dynamical systems involving bosons and fermions (for a 
review, see Ref. Q). They play a role in diverse phenom- 
ena including nuclear and molecular spectroscopy [Tol — 
[l2j . quantum phase transitions [l3T ] and mixed regular 
and chaotic dynamics [13]. Here we consider the SU(3) 
symmetry in view of its significance for deformed nuclei, 
as reco gniz ed in the Elliott and symplectic shell mod- 
els [15|, ll6( and the IBM. We use the mathematical al- 
gorithm to construct, order b y o rder, all possible inter- 
actions with a given PDS [H1II1, apply it to the SU(3) 
limit of the IBM, and illustrate with a concrete example 
how the PDS and data constrain the form and strength 
of higher-order interactions. 

The IBM describes low-energy collective states of the 
nucleus in terms of N monopole (s) and quadrupole (d) 
bosons representing pairs of nucleons. The dynamical al- 
gebra is U(6) with generators in terms of which operators 
of all physical observables can be written. The classifica- 
tion of states in the SU(3) limit is [19( 



U(6) D SU(3) D SO(3) D SO(2) 
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where underneath each algebra the associated labels are 
given (K is a multiplicity label needed in the SU(3) D 
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SO(3) reduction). These define the Elliott basis 
\[N](\, /i)KLM) , from which the Vergados basis 
\[N](X, fi)xLM) , is obtained by a standard orthogonal- 
ization procedure. The classification (JTJ) assumes a sym- 
metric U(6) irreducible representation (irrep) [N] which 
is appropriate for the IBM. Apart from terms involv- 
ing the conserved total boson number operator N, a 
rotational-invariant Hamiltonian with SU(3) DS has the 
form 



Hi 



DS 



aiC 2 [SU(3)]+a 2 C 2 [SO(3)]+a 3 C 3 [SU(3)] , (2) 



where C n [G] is the n order Casimir operator of the Lie 
algebra G and a% are coefficients. This form exhausts 
all independent Casimir operators of SU(3) and SO(3), 
that is, any other commuting operator can be written as a 
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TABLE I: Normalized two- and three-boson SU(3) tensors. 
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function of those appearing in Eq. ([2]) . Hds is completely 
solvable with eigenenergies 

Ev S =a 1 f 2 (\, f j,)+a 2 L(L + l) + a :i M\, fJ ,) , (3) 

where / 2 (A,/i) = A 2 + (A + /x)(/x + 3) and / 3 (A,/x) = 
(A — /i)(2A + /i + 3)(A + 2/i + 3). The spectrum resembles 
that of a quadrupole axially-deformed rotor with eigen- 
states arranged in SU(3) multiplets and K corresponds 
geometrically to the projection of the angular momen- 
tum on the symmetry axis. The Hamiltonian -Hds is 
genuinely many-body (with interactions that are up to 
third order in the bosons). Its applicability is limited, 
however, since only three independent operators exist, 
and states in different -ftT-bands with the same (A,/z)L 
are degenerate. Flexibility can be considerably increased 
by introducing interactions with PDS. The method to 
construct such interactions is based on an expansion of 
the Hamiltonian in terms of tensors which annihilate pre- 
scribed set of states [13, EH]. In the present study, the 
tensors involve n-boson creation and annihilation opera- 
tors with definite character under the SU(3) chain ([TJ, 

B[n](\,p,)xim> An 5 ]0,A)xfm = ( _ )"' (^[n] ( \,fi)xt ,-m) t ' W 

The SU(3) tensor operators for n=2 and 3 are given in 
Table HI Of particular interest are the operators with 
(A, fi) ^ (2n, 0) because the corresponding annihilation 
operators yield zero when acting on the ground-band 
members \[N](2N,0)K = 0,LM) (and possibly other 
states). Interactions involving these operators can be 
added to the Hamiltonian without destroying solv- 
ability of part of its spectrum. Two such operators, Pq 
and P 2m , exist for n=2, and allow the construction of 



TABLE II: Number of interactions in the IBM. 
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"On the left of >-> is the number of interactions of a given order; 
this reduces to the number on the right of i-> if one is only interested 
in excitation energies in a single nucleus. 



an IBM Hamiltonian with up to two-boson interactions 
that have a solvable ground band \[N](2N,0)K = 0, LM) 
and a solvable 7 band |[iV](2iV - 4,2)K = 2, LM). A 
two-body Hamiltonian with SU(3) PDS can be applied 
to 168 Er and the excellent SU(3) description of the en- 
ergies and E2 properties of these bands can be retained 
while lifting the degeneracy of the /3 and 7 bands 10 . 

The most general (2 + 3)-body Hamiltonian with SU(3) 
PDS can be written in terms of the operators given in 
Table H 

#pds = h 2 P^p2 + h P^P +g 4 WlW i + g 3 WlW 3 
+92V2 ■ V2 + g\Wl ■ W 2 + g c 2 (Vl ■ W 2 + W\ ■ V 2 ) 
+ 3o a AtA + gfrvtWo + g c (A f W + W$A), (5) 

with 2 + 8 interactions strengths hi and gf . Terms in- 
volving the operator C , 2[SO(3)] = L 2 can be added to 
this Hamiltonian, as is done in Eq. (|2|). To illustrate the 
increase in flexibility of a Hamiltonian with SU(3) PDS, 
we list in Table |TT] the number of interactions under the 
different scenarios. Up to third order, a general rotation- 



3 



> 



OO 

c 1 



- 14 + 


10 + 10 + 

9 + 


156^ , 

64 Lrd 92 

14 + io + 

10 + 9 + 


— w — 10+ 

10 + 9 + 

8 + 


19 + 

10 + 


8 + 8 + 

— T 

— 6 + — §: 

4+ _4+ 

^— 0+ 2 + 


12 + 8 + 

— 10 _ 4: —4: 

2 + 4+ 

— 8 / r 

6 + 


17+ 8 + 7+ " 

— 10 + _ 4 + —it 
2 + 4+ 

A+ Z 

— 8+ ° 


8 + ' 

6 + 


(2N-4.2) 


(2N-4.2) 
6 + 


4 t 

2 + 




4 + 

2+ 


4 t 


- — + 
(2N,0) 


SU(3)DS 


~° Experiment 


( =|* SU(3)PDS- 



FIG. 1: Observed spectrum of 156 Gd |25[| compared with the calculated spectra of Hx>s © with SU(3) DS and of -ffp DS © 
with SU(3) PDS for N = 12 and parameters given in the text. Bands are labeled with the SU(3) quantum numbers (A,/i). 



ally invariant Hamiltonian has 26 independent interac- 
tions, decreasing to 16 if one is only interested in exci- 
tation energies in a single nucleus. (This excludes terms 
involving N). A Hamiltonian with SU(3) DS has, up to 
third order, 8 independent terms but 5 of them (TV, N 2 , 
N 3 , NL 2 , and iVC^ [SU(3)]) are constant in a single nu- 
cleus or can be absorbed in an interaction of lower order, 
leaving only the 3 genuinely independent terms shown in 
Eq. ([5]). The corresponding numbers for a Hamiltonian 
with SU(3) PDS are 15 and 9. The latter number agrees 
with the 10 terms in the Hamiltonian ([5]) which lacks L 2 
but includes the combinations NPl-P 2 and NPqPq. We 
conclude from Table HI1 that more than half of all possible 
interactions in the IBM have in fact an SU(3) PDS. 

Several SU(3)-preserving interactions are contained in 
the expression ©. Specifically, 9 2 = 2N(2N + 3) - C 2 
corresponds to h = h 2 = !8; (N-2)9 2 : g$ = 54, g„=g 2 = 
g%=g 3 = , 9 4 = 30; C 3 + (2JV+3) [3(9 2 -2 A(4 A+3)] : 5 g = 648 
and Cl - (AN + 3)L 2 : h 2 = -108, g% = 9 5 g = — 3c?g = 270, 
.92/5 = S2/ 21 = 52/\/l05 = 24/13, 54 = -120. The three 
terms involving C*„[SU(3)] are included in H^s ©■ The 
(integrity basis) term Cl — —4\^3Q-(L x L)( 2 ' is com- 
posed of SU(3) generators, hence is diagonal in (A,/x), 
but breaks the if-degeneracy of the exact DS. Its impact 
on nuclear spectroscopy has been well studied in the sym- 
plectic shell model and the IBM HH-il. The PDS notion 
goes a step further by allowing SU(3) mixing in most (but 
not all) of the eigenstates of the Hamiltonian. 

As noted, in general, Hpvs © does not preserve 
SU(3) yet, by construction, for any choice of parame- 
ters the ground-band members \[N](2N,0)K = 0,LM) 
are solvable. For specific choices, additional solvable 
states are obtained. In particular, by choosing only the 
^0, 9o, So) 9o terms and H DS ©, the states \[N](2N - 



A,2k)K = 2k, LM), k = 1, 2, . . .(among which the 7- 
band members with fc = l) remain solvable with energies 
-Eds © • This case therefore has the same solvable states 
as the two-body Hamiltonian with SU(3) PDS consid- 
ered in Ref. [l(| ; the additional three-body terms lead to 
a different mixing of the non-solvable states. 

Another class of Hamiltonians with SU(3) PDS exists 
which has solvable /3-band members |[iV](2jV — 4, 2)K = 
0, LM) with energies -Eds ©• This follows from the 
structure of the relevant Hamiltonian, 

#pds = #ds + mW\-W 2 + n 3 wlw 3 , (6) 
and the fact that W 2m and W 3m annihilate the in- 
trinsic state of the band, |/3) cx (V2Pq - P2o)( st + 
\Z2do) N ~ 2 \0}. The property of solvability of the ground 
and /3 bands can be exploited in the following way. The 
Hamiltonian Hos © has a rotor spectrum with char- 
acteristic L(L + 1) splitting for all bands. Deviations 
from this pattern are often observed for the 7 band of 
deformed nuclei and are indicative of 7-soft or triaxial 
behavior [24[. We illustrate the procedure with an ap- 
plication to 156 Gd. The parameters ot\ = —7.6 keV, 
a 2 = 12.0 keV, and a 3 = in Hus are fixed from the 
excitation energy of the /3-band head and the moments 
of inertia of the ground and fi bands. This completely 
determines the SU(3) DS spectrum, shown on the left 
of Fig. [TJ which is characterized by degenerate j3 and 7 
bands. In the observed spectrum these bands are not 
degenerate and, more importantly, the 7-band energies 
display an odd-even staggering. This effect can be visu- 
alized by plotting the quantity (26j 



Y{L)= 2_L-l x E(L)-E(L-i; 



1, 



(7) 



L " E(L) - E(L — 2 / 

where E(L) is the excitation energy of a 7-band level 
with angular momentum L. For a rotor this quantity is 
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FIG. 2: Observed and calculated (SU(3) DS and PDS) odd- 
even staggering of the 7 band in 156 Gd. 
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flat, Y(L) = 0, as illustrated in Fig. H with the SU(3) DS 
calculation. The data, however, show considerable odd- 
even staggering which can be well described by a combi- 
nation of three-body interactions with 772 = —18.1 keV 
and 773 = 46.2 keV. The calculated staggering increases 
with L which agrees with the experiment up to L = 10. 
For L > 10 the observed staggering changes character, 
a phenomena requiring higher angular momentum pairs, 
which are beyond the scope of the standard (s, d) IBM 
description. The two interactions W\ ■ W2 and W$ ■ W3 
induce a mixing of the 7 band with higher-lying excited 
bands. Other approaches advocating the coupling of the 
7 band to the (3 band [23[ or to the ground band [l?} 
fail to describe the odd-even staggering in 156 Gd. For 
the PDS calculation, the wave functions of the states in 
the 7 band involve 15% SU(3) admixtures into the domi- 
nant (27V — 4, 2) component. Higher bands exhibit larger 
SU(3) mixing and their wave functions are spread over 
many SU(3) irreps, as shown for the K = O3 band in 
Fig. 3. This complex SU(3) decomposition is in marked 
contrast to the SU(3)-purity of the ground (K = 0i) and 
P (K = O2) bands. Such strong symmetry-breaking can- 
not be treated in perturbation theory. 

Previous studies of triaxiality in the IBM framework 
have employed only the cubic 773 term of Eq. © [28|, [29j ■ 
The current work hints that both 772 and 773 terms are 
necessary for an accurate description of odd-even stag- 
gering in deformed nuclei. This highlights the capacity 
of the PDS approach to identify novel relevant terms of a 
given order. We emphasize that the PDS results for the 
7 band are obtained without altering the good agreement 
for the ground and /3 bands, already achieved with the 
SU(3) DS calculation. This is further illustrated with the 
E2 transitions in 156 Gd. The observed B(E2) values be- 
tween ground, f3, and 7 bands are shown in Table Hill and 
compared to the results of the SU(3) DS and PDS calcu- 
lations. The effective boson charge eb = 0.166 eb in the 
electric quadrupole operator eh[s^d+d^s + x(d^d)^] and 
the value \ = -0.168 arc fitted to the S(E2; 2+ 0+) 
and -B(E2; 2 + — > 4+) values. The E2 transitions between 
ground and f3 bands can be calculated analytically (30j . 



FIG. 3: SU(3) decomposition of wave functions of L — 
states in the K = 0i, O2, O3 bands for the PDS calculation. 



TABLE III: Observed and calculated B(E2) values in 156 Gd. 
The parameters of the E2 operator are given in the text. 
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and remain valid in SU(3) PDS. Transitions involving 7- 
band members are different in SU(3) DS and PDS, and 
are computed numerically for the latter. It is seen from 
Table Hill that the mixing of the 7 band with higher- lying 
excited bands improves the agreement with the data in 
most cases. 

In summary, we have identified several classes of (2+3)- 
body IBM Hamiltonians with SU(3) PDS, and obtained 
an improved description of signature splitting in the 7 
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band of 156 Gd. The analysis serves to highlight the mer- 
its gained by using the notion of PDS as a tool for se- 
lecting higher-order terms in systems where a prescribed 
symmetry is not obeyed uniformly. On one hand, the 
PDS approach allows more flexibility by relaxing the con- 
straints of an exact DS. On the other hand, the PDS picks 
particular symmetry-breaking terms which do not de- 
stroy results previously obtained with a DS for a segment 
of the spectrum. The PDS construction is implemented 
order by order, yet the scheme is non-perturbative in 
the sense that the non-solvable states experience strong 
symmetry-breaking. These virtues can be exploited in 
attempts to extend the ab-initio and beyond-mean-field 



methods to heavy nuclei. The present work motivates 
and sets the stage for further exploring the impact of 
PDS with higher-order terms on the dynamics in quan- 
tum many-body systems. 
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